NASA/TM -^7. 


208146 


V. K. Garg 

AYT Corporation. 


R. E. Gaugler 

Turbine Branch, 
Turbomachinery & Propulsion 
Systems Division. 

NASA Lewis Research Center, 
Mail Stop 5-11, 
Cleveland, OH 44135 


At/S 


Effect of Velocity and 
Temperature Distribution 
Hole Exit on Film Cooling 
Turbine Blades 

An existing three-dimensional Navier-Stokes code (Amone et al, 1991), modified 
to include film cooling considerations (Garg and Gaugler, 1994), has been used to 
study the effect of coolant velocity and temperature distribution at the hole exit on 
the heat transfer coefficient on three film-cooled turbine blades, namely, the C3X 
vane, the VK1 rotor , and the ACE rotor . Results are also compared with the experi- 
mental data for all the blades. Moreover, Mayle’s transition criterion (1991), For- 
est's model for augmentation of leading edge heat transfer due to free -stream turbu- 
lence ( 1977), and Crawford's model for augmentation of eddy viscosity due to film 
cooling (Crawford et al, 1980) are used. Use of Mayle's and Forest's models is 
relevant only for the ACE rotor due to the absence of showerhead cooling on this 
rotor. It is found that, in some cases, the effect of distribution of coolant velocity and 
temperature a t th e hole exit can be as much as 60 percent on the heat transfer 
coefficient at the blade suction surface, and 50 percent at the pressure surface. Also, 
different effects are observed on the pressure and suction surface depending upon 
the blade as well as upon the hole shape, conical or cylindrical. 
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1 Introduction 

There is a growing tendency these days to use higher and 
higher temperatures at the inlet to a turbine since it yields higher 
thermal efficiency. Modem gas turbine engines are designed to 
operate at inlet temperatures of 1800-2000 K, which are far 
beyond the allowable metal temperatures. Thus, to maintain 
acceptable life and safety standards, the structural elements need 
to be protected against the severe thermal loads. This calls for 
an efficient cooling system. One such cooling technique cur- 
rently used for high-temperature turbines is film cooling. In 
this technique, cooler air is injected into the high-temperature 
boundary layer on the blade surface. Since the cooler air is 
bled directly from the compressor before it passes through the 
combustion chamber, it represents a loss in the total power 
output. The designer’s goal is therefore to minimize the coolant 
necessary to insure adequate turbine life. 

To this end, considerable effort has been devoted to under- 
standing the coolant film behavior and its interaction with the 
mainstream flow. The film cooling performance is influenced 
by the wall curvature, three-dimensional external flow structure, 
free-stream turbulence, compressibility, flow unsteadiness, the 
hole size, shape and location, and the angle of injection. Many 
studies on film cooling have been confined to simple geometries, 
for example, two-dimensional flat and curved plates in steady, 
incompressible flow. A survey of work up to 1971 has been 
provided by Goldstein (1971). While several further studies in 
this field have been summarized by Garg and Gaugler ( 1993, 
1994, 1996), we will discuss some relevant ones here from the 
point of view of the present study. 

Bergeles et al. (1980) devised a finite-difference code with 
a semi-elliptic treatment of the flow field in the neighborhood 
of the injection holes. The comparison of predicted results with 
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experimental data on a flat plate showed poor agreement in the 
vicinity of the hole partly due to the assumption of a uniform 
coolant velocity at the hole exit. Schdnung and Rodi (1987) 
presented a two-dimensional boundary layer model with a modi- 
fication for three-dimensional elliptic flow for simulating the 
effects of film cooling by a single row of holes. The region near 
the injection was leapt over and new profiles were prescribed at 
a certain distance behind the blowing region. Later, Haas et al. 
( 1992) extended Schonung and Rodi’s ( 1987) model to account 
for density differences between the hot gas and the injected 
coolant gas. However, both the models did not account for 
the effects of curvature and multiple rows of holes. Tafti and 
Yavuzkurt (1990) developed a two-dimensional injection 
model for use with a two-dimensional low-Reynolds-number 
k - e model boundary layer code for film-cooling applications. 
They specified a uniform coolant temperature but nonuniform 
velocity profile at the hole exit. Dibelius et al. ( 1990) developed 
an elliptic procedure near the injection area but a partially para- 
bolic procedure far downstream for film cooling on a flat plate. 
They prescribed a uniform coolant velocity but slightly nonuni- 
form coolant temperature at the hole exit. 

Vogel (1991) described a three-dimensional Navier-Stokes 
code for film-cooled gas turbine blades that couples the flow 
problem over the blade with the three-dimensional heat conduc- 
tion problem within the blade. Using a uniform coolant velocity 
at the hole exit, Vogel compared the heat transfer coefficients 
with experimental data for injection on a flat plate. He could 
not get reliable results on a turbine blade due to the computer 
storage limitations that did not allow more than four to six mesh 
cells within the hole opening on the blade surface. Benz and 
Wittig (1992) analyzed the elliptic interaction of film-cooling 
air with the main flow by simultaneously computing the coolant 
and main flows for film cooling at the leading edge of a turbine 
blade. They prescribed a uniform coolant flow at entry to the 
injection pipe, and found a nonuniform velocity profile at the 
hole exit. They did not provide any information about the tem- 
perature field. Amer et al. ( 1992) compared two forms each of 
the k- e and k- uj family of turbulence models for film cooling, 
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and found all of them to be inappropriate. Realizing the impor- 
tance of coolant distribution at the hole exit, they tried a linear, 
a parabolic, and the ^th power-law velocity profile for the cool- 
ant at the hole exit. Finding only small differences due to the 
parabolic and power-law profiles, they provided most results 
using the latter profile. 

Domey and Davis (1993) analyzed the film-cooling effec- 
tiveness from one and two rows of holes on a turbine vane, 
using Rai’s ( 1989) numerical technique. They carried out both 
two- and three-dimensional simulations, but represented each 
hole by just two grid points, and specified uniform velocity 
and temperature of the coolant at the hole exit. Choi (1993) 
developed a multizone code for film cooling with a fine mesh 
near the coolant hole and injection pipe overlapped on the global 
coarse grid. He specified uniform conditions for the coolant at 
entry to the injection pipe but provided no details for the coolant 
profiles at the hole exit. Leylek and Zerkle ( 1994) analyzed the 
coupled problem of flow in the plenum, film-hole, and cross- 
stream region for film cooling on a flat plate using the k-c 
model of turbulence in a three-dimensional Navier- Stokes 
code. They found the flow within the film-hole to be extremely 
complex, containing counterrotating vortices and local jetting 
effects that make the flow field in this region highly elliptic. 
The distribution of dependent variables at the hole exit plane 
resulted from the interaction of three competing mechanisms, 
namely, counterrotating structure and local jetting effects within 
the film-hole, and crossflow blockage. 

Fougeres and Heider (1994) solved the unsteady three-di- 
mensional Navier- Stokes equations, completed by a mixing- 
length turbulence model, using a finite volume technique. They 
presented two applications of the multidomain code: one for a 
single row of hot jets injected into a flat-plate turbulent bound- 
ary layer, and another for a plane nozzle guide vane with two 
rows of staggered holes on the pressure as well as suction side 
of the vane. They specified uniform conditions for the coolant 
at entry to the injection pipe but provided no details for the 
coolant profiles at the hole exit. Weigand and Harasgama 
( 1994) carried out a numerical investigation of film cooling on 
a turbine rotor blade using the Dawes ( 1993) code that utilizes 
an unstructured solution adaptive grid methodology for solving 
three-dimensional Navier- Stokes equations. They discretized 
the plenum chamber and the injection pipe. However, a rather 
academic case of blowing in the tangential direction was studied 
due to limitations of the code. As such, comparison with experi- 
mental data was not possible. 

Hall et al. ( 1994) analyzed the showerhead film cooling on 
the C3X vane with a multiblock, three-dimensional Navier- 


Stokes code using the Baldwin -Lomax turbulence model. Tak- 
ing advantage of the spanwise periodicity of the planar C3X 
vane, the computational span was restricted to just one spanwise 
pitch of the showerhead holes. The final airfoil C-grid had over 
2.1 million grid points with a 17 x 17 grid patch on each of 
the five holes. Such a grid refinement cannot be handled by 
present-day computers if the vane were annular, requiring the 
whole span to be analyzed, not just a slice of it. The computa- 
tional domain included the injection pipe but not the plenum 
chamber. They specified uniform coolant conditions at entry to 
the injection pipe but provided no details of the coolant distribu- 
tion at the hole exit. 

This survey indicates that the effect of coolant distribution 
at the hole exit, though important, has not yet been investigated 
on real turbine blades. This is precisely the objective of the 
present study. Herein, an existing three-dimensional Navier - 
Stokes code (Amone et al., 1991), modified to include film- 
cooling considerations (Garg and Gaugler, 1994), has been 
used to study the effect of coolant velocity and temperature 
profiles at the hole exit on the heat transfer characteristics of 
film-cooled turbine blades. Comparison with experimental data 
for a C3X vane with nine rows of film cooling holes (Hylton 
et al., 1988), for a VKI rotor with six rows of holes (Camci 
and Arts, 1990), and for the ACE rotor with three rows of holes 
(Norton et al., 1990), is provided. 

2 Analysis 

The three-dimensional Navier- Stokes code of Amone et al. 
( 1991 ) for the analysis of turbomachinery flows was modified 
by Garg and Gaugler (1994) to include film-cooling effects 
and Mayle’s transition criterion (Mayle, 1991). It was further 
modified to include Forest’s model for augmentation of leading 
edge heat transfer due to free-stream turbulence (Forest, 1977), 
and Crawford’s model for augmentation of eddy viscosity due 
to film cooling (Crawford et al., 1980). Briefly, the code is an 
explicit, multigrid, cell-centered, finite volume code with an 
algebraic turbulence model. The Navier- Stokes equations in a 
rotating Cartesian coordinate system are mapped onto a general 
body-fitted coordinate system using standard techniques. Vis- 
cous effects in the streamwise direction are neglected in compar- 
ison to those in the other two directions. Justification for this 
assumption is provided by Garg and Gaugler (1994). The 
multistage Runge-Kutta scheme developed by Jameson et al. 
(1981) is used to advance the flow solution in time from an 
initial guess to the steady state. A spatially varying time step 
along with a CFL number of 5 was used to speed convergence 


Nomenclature 


a y b = semi-axes of the ellipse describ- 
ing the hole exit on the blade 
surface 

<a 0 -a 3 = coefficients in Eq. ( 1 ) 

B f = blowing ratio - (p c V' c )/(p 3O ?7 0C ) 
c = true chord of the blade 
C F = coefficient in Eq. (2) 
d - coolant hole diameter 
h = heat transfer coefficient based 
on ( T 0 - T w ) 

h a ~ standard value = 1135.6 W/m 2 - 
K = 200 Btu/hr-ft 2 -R 
k = von Karman constant 
l = mixing length 
L = length of the injection pipe 
M = Mach number 
p = pressure 

PD = penetration distance (Craw- 
ford’s model) 


r = coolant hole radius = d/2 
Re = Reynolds number based on the 
true chord length 

s = distance from the leading edge 
along the pressure or suction sur- 
face 

s 0 , Zo = coordinates of the center of hole 
exit on the blade surface 
T = temperature 
Tu = turbulent intensity 
U = main-flow velocity 
V = coolant velocity at any point 
within the hole exit 
x = streamwise distance from the 
point of injection 

y + = dimensionless distance of the 
first point off the blade surface 
z = distance along the span 
6 = boundary layer thickness 


= eddy diffusivity for momentum 
y = ratio of specific heats 
77 , £ = coordinates from the center of 
hole exit (see Eq. (1)) 

6 = momentum thickness 
1 / = kinematic viscosity 
p = density 

Subscripts 

1 = at inlet 

2 = at exit 

c — for coolant (average value) 
e = edge of boundary layer 
m = maximum value 
n = corresponding to uncooled blade 
o = stagnation value 
w = at the blade surface 
oc = local free-stream value 
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to the steady state. Eigenvalue -scaled artificial dissipation and 
variable-coefficient implicit residual smoothing are used along 
with a full-multigrid method. 

The effects of film cooling have been incorporated into the 
code in the form of appropriate boundary conditions at the hole 
locations on the blade surface. Each hole exit (generally an 
ellipse with semi-axes a and b in the plane of the blade surface) 
is represented by several control volumes (about 20) having a 
total area equal to the area of the hole exit, and passing the 
same coolant mass flow. Different velocity and temperature 
profiles for the injected gas can be specified at the hole exit. 
Earlier studies indicate that typical gas turbine film cooling 
situations involve steep polynomial profiles with their maxima 
shifted in upstream or downstream direction from the jet center- 
line. Based on a detailed computation through the plenum cham- 
ber and injection pipe, Leylek and Zerkle ( 1994) found that for 
high Lid (>3.0) and high blowing ratio (^1.0) or for low L! 
d and low blowing ratio, the exit velocity profile is very much 
like a fully developed turbulent pipe flow profile. For other 
combinations of Lid and blowing ratio, Leylek and Zerkle 
(1994) found that the exit velocity profile is polynomial with 
shifted maxima. For the cases reported here, turbulent (|th 
power-law) and polynomial profiles of the form 

y = O} Y _ T “ = (1 - v 2 - C 2 )(a° + a, T) + a 2 ri 2 ) (1) 
where 


7? « (s - s 0 )la , C = (Z “ Zo)fb 

were specified, in order to study their effect on the heat transfer 
coefficient at the blade surface. Note that Eq. ( 1 ) automatically 
satisfies V = 0 condition on the boundary of the hole exit, and 
in it ( s Q , Zo ) denotes the center of hole exit, and rj and £ are 
interchangeable. Using the coolant mass conservation condition, 
coefficients a oy a x , a 2 can be related to V m , the maximum value 
of V and its location, V,. Then using the coolant energy conser- 
vation, (3 3 can be related to the other coefficients. These relations 
are 


a 2 = 


1 

1 + 6 Vj 


12 - 


6V„(1 - 3V, 2 ) ' 
(1-V ?) 2 . 


2V,V„ „ „ „ 

a ' = (i - vjy - 2a2Vl ’ a ° = 2 ~ a * /6 ’ 

a 3 = all2 + ( a] + 20^2)^24 + 02 / 8 O 

For all three blades analyzed here, the experimentally deter- 
mined temperatures were specified at the blade surface, and wall 
heat flux was calculated. The algebraic mixing length turbulence 
model of Baldwin and Lomax (1978) was used. This model 
has been used satisfactorily by Boyle and Giel (1992), Ameri 
and Amone (1994, 1996), and Boyle and Ameri (1997) for 
heat transfer calculations on turbine blades without film cooling, 
and by Hall et al. (1994), and Garg and Gaugier (1994) with 
film cooling. In fact, Ameri and Amone (1994) compared the 
Baldwin-Lomax model and Coakley’s q-u model against ex- 
perimental data of Graziani et al. (1980), and found that the 
algebraic model was able to produce many of the flow features 
better than the two-equation model. They further state that this 
conclusion is strengthened when one takes into account the 
relative economy of computations with the algebraic model. 

Moreover, Mayle’s transition criterion (1991) has been im- 
plemented in the code. The incoming flow in the experimental 
tests on the ACE rotor (Norton et al., 1990) had a turbulence 
intensity of 4.0 percent, and it is assumed constant for applica- 
tion of Mayle’s model. For the C3X vane and the VKI rotor, 
Mayle’s criterion is immaterial since the flow turns turbulent 
at the leading edge itself due to showerhead injection. 


2.1 Forest’s Model. The strong favorable free-stream 
pressure gradients near the leading edge result in laminar flow, 
even with high free-stream turbulence. The laminar heat transfer 
at the leading edge is, however, increased by the free-stream 
turbulence. The model of Forest (1977) was incorporated in 
order to account for this increase. According to this model 

e M = CjdTuJJ^ (2) 

where the mixing length, /, is defined as 

/ = min (ky, 0.08<5) 

While several ways to calculate the boundary layer thickness 
were tried following Davis et al. ( 1988), it was found that the 
best method was to relate it to the momentum thickness whose 
calculation is described later on. The coefficient C F is found 
from 

C 0-75 

F 1 + 0.04/ 1 X | 

where 

v ds 


The momentum thickness is calculated using Thwaites’ method, 
as given by White (1974). According to this, the momentum 
thickness at a location s 1 along the blade surface is given by 


0 2 (s ,) 


0A5v r> 

UtM Jo 


u]ds 


The origin in this calculation is taken as the geometric stagna- 
tion point. This calculation of e M due to free-stream turbulence 
was carried out only when the flow was laminar. Otherwise, 
the analysis would give significantly increased heat transfer for 
fully turbulent flow (Boyle, 1991 ). For the cases analyzed here, 
it implies that Forest’s model is used only for the ACE rotor 
since for both the VKI rotor and the C3X vane, the flow is 
turbulent at the leading edge itself due to the shower-head injec- 
tion. 


2.2 Crawford’s Model. Crawford’s model for augmen- 
tation of eddy viscosity due to film cooling (Crawford et al., 
1980) was developed for use in a two-dimensional boundary 
layer code to predict heat transfer in full-coverage film cooling 
situations on flat surfaces. Though not really applicable to the 
situation, it was used by Stepka and Gaugier (1983 ) to compare 
the predicted Stanton numbers with the experimental data for 
rows of film-cooling holes around the circumference of a cylin- 
der in crossflow. In the absence of a better model, we adopted 
it also for the present calculations. According to this model, the 
eddy diffusivity for momentum in the presence of film cooling is 
modeled by algebraically augmenting the Prandtl mixing length 
using 



where the “3-D” subscript refers to the three-dimensional mix- 
ing length from the Baldwin and Lomax model, and the “0” 
denotes the departure due to jet-boundary layer interaction, 
given by 


F-f, 


( 0 - 
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/= exp[-(x/<5)/2], 

Knax* = 0.0353 exp(2.65B r ) for slant angled injection 

\nux^i = 0.0177 exp(2.64B r ) for compound angled injection 

Km* = 0.0601 exp (3. 46 B r ) for normal injection 

where x is the streamwise distance measured from the point of 
injection, and PD is the penetration distance. While Crawford’s 
model is essentially two dimensional, we extended it to three 
dimensions by linearly augmenting the eddy viscosity in the 
spanwise direction between holes. 

3 Blades and Experimental Details 

The experimental data on the C3X vane have been provided 
by Hylton et al. (1988). Figure 1 shows the C3X vane with 
the film-cooling hole details. The test vane was internally cooled 
by an array of ten radial cooling holes (not shown in Fig. 1) 
in the active part of the vane. No heat transfer measurements 
were made in the actual film-cooled nose piece of the vane as 
it was thermally isolated from the rest of the vane. In addition 
to the details available in Fig. 1, the two rows of holes on the 
pressure surface were located at s/s m — 0.197, 0.225, and on 
the suction surface at s/s m = 0.254, 0.276. One row of show- 
erhead holes was located at the geometric leading edge, with 
one row on the suction side and three on the pressure side. The 
active part of the test vane surface was instrumented with 123 
0.51-mm-dia sheathed CA thermocouples to measure the tem- 
perature at the blade surface. The thermocouple junctions were 
located in a plane 2.54 mm off midspan. The heat transfer 
coefficient for each radial cooling hole was calculated from the 
hole diameter, measured coolant flow rate, and coolant tempera- 
ture with a correction applied for thermal entry region effects. 
The internal temperature field of the test vane was obtained 
from a finite element solution of the steady-state heat conduction 
equation, using the measured surface temperatures as boundary 
conditions. Experimental values of the heat flux at the blade 
surface were then obtained from the normal temperature gradi- 
ent at the blade surface, and these were used to derive the local 
heat transfer coefficient at a spanwise location 2.54 mm off the 
midspan of the blade. More details are available from Hylton 
etal. (1988). 

The experimental data on the VKI rotor have been provided 
by Camci and Arts ( 1990), using the short-duration VKI Isen- 
tropic Compression Tube facility. Figure 2 shows the VKI rotor 
geometry along with cooling hole details. Three staggered rows 
of cylindrical cooling holes (d = 0.8 mm; s/c = -0.0285, 0, 
0.0285) were located around the leading edge. The row and 
hole spacings were both 2.48 mm. These holes were spanwise 
angled at 30 deg from the tangential direction and drilled in a 


2 ROWS OF HOLES 
35* ANGLE CHORDWISE 
3d SPANWISE SPACING 



$ ROWS OF STAGGERED HOLES 
45* ANGLE SPANWISE 
7Sd SPANWISE SPACING 


HOLE DIA. d « 0.99 mm 
4d ROW SPACING 

HOLES NORMAL TO DIRECTION UNSPECIFIED 
Fig. 1 C3X vane and cooling hole details 



Fig. 2 VKI rotor and cooling hole details 


plane perpendicular to the blade surface. Two staggered rows 
of conical holes (d = 0.8 mm; s/c - 0.206, 0.237) were located 
on the suction side. Hie row and hole spacings were respectively 
2.48 and 2.64 mm. These holes were inclined at 37 and 43 deg 
with respect to the local blade surface and drilled in a plane 
perpendicular to the span. One row of conical holes (d = 0.8 
mm; s/c = —0.315) was located along the pressure side. The 
hole spacing was 2.64 mm. These holes were inclined at 35 
deg with respect to the local blade surface and drilled in a plane 
perpendicular to the span. The blade instrumented for heat flux 
measurements was milled from 4 ‘Macor’ * glass ceramic and 45 
platinum thin films were applied on its surface. Three indepen- 
dent cavities were drilled along the blade height to act as plenum 
chambers. The local wall convective heat flux was deduced 
from the corresponding time-dependent surface temperature 
evolution, provided by the platinum thin-film gages. The wall 
temperature/wall heat flux conversion was obtained from an 
electrical analogy, simulating a one-dimensional semi-infinite 
body configuration. More details are available from Camci and 
Arts (1990). 

The experimental data on the ACE rotor have been provided 
by Norton et al. ( 1990), using the Isentropic Light Piston Tun- 
nel at Oxford University. The blade profile tested was the mid- 
height streamline section of a highly loaded, transonic, research 
turbine rotor. Build X of the ACE rotor, shown in Fig. 3, was 
used for comparison purposes. In this, one row of cylindrical 
holes (d = 0.65 mm; s/s m = 0.7) was located on the suction 
surface, with a hole spacing of Ad, Two staggered rows of 
cylindrical holes ( d = 0.65 mm; s/s m = 0.5, 0.56) were located 
on the pressure surface, with a row and hole spacing of Ad 
each. All holes were inclined 30 deg to the local blade surface 
and drilled in a plane perpendicular to the span. The instru- 
mented blades were manufactured from “Macor” glass ceramic 
onto which were placed thin film platinum resistance thermome- 
ters and pads. The heat transfer measurement technique is simi- 
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lar to that used by Camci and Arts (1990). More details are 
available from Norton et al. ( 1990). 

4 Computational Details 

For all the blades, the computational span, shown in Fig. 4, 
is different from the real span. The ordinate in each part of Fig. 
4 denotes the distance along the blade surface in the spanwise 
direction, while the abscissa denotes the distance along the blade 
surface in the chordwise direction, both normalized by the hole 





Fig. 4 Computational span for the C3X vane, VKI rotor, and the ACE 
rotor 


Table 1 Parameter values for the cases analyzed 
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Film Cooling Parameters (Derived) 
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radius, r. For the VKI rotor, the hole spacing in the span direc- 
tion was assumed to be 2.56 mm (=6.4r) for all holes, repre- 
senting a deviation of ±3 percent from the actual values of 2.48 
and 2.64 mm. It may be noted that the abscissa in Fig. 4 has 
breaks so as to accommodate all the rows of holes. The shape 
and orientation of the hole openings in Fig. 4 are a direct conse- 
quence of the angles the holes make with the spanwise or 
chordwise direction. For all the blades, the pattern of holes 
shown in Fig. 4 is repeated in the spanwise direction. For the 
C3X vane and the VKI rotor, periodic boundary conditions are 
imposed due to showerhead injection, while for the ACE rotor, 
symmetry conditions are imposed on the ends of the computa- 
tional span. 

Since the hole diameter on all the blades is less than 1 mm, 
the grid size has to be varied along the blade chord. For compu- 
tational accuracy, the ratio of two adjacent grid sizes in any 
direction was kept within 0.76 to 1.3. A periodic C-grid with 
up to one million grid points was used. For the C3X vane, the 
grid used was 281 X 45 X 81, while for the VKI and ACE 
rotors, it was 353 X 53 X 17 and 313 X 53 x 17, respectively, 
where the first number represents the number of grid points 
along the main flow direction, the second in the blade-to-blade 
direction, and the third in the span direction. Normal to the 
blade surface is the dense viscous grid, with y * < 1 for the first 



Fig. 5 Coolant velocity and temperature profiles at the hole exit 
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Fig. 6 Blade surface temperature for various cases (C3X vane) Fig. 8 Effect of coolant velocity and temperature profiles on the normal- 

ized heat transfer coefficient at the C3X vane surface (case 44155) 


point off the blade surface, following Boyle and Giel (1992). 
Computations were run on the 8-processor Cray Y-MP super- 
computer at the NASA Lewis Research Center, and on the 16- 
processor C-90 supercomputer at NASA Ames Research Cen- 
ter. The code requires about 60 million words (Mw) of storage 
and takes about 20 seconds per iteration (full-multigrid) on the 
C-90 machine for one million grid points. For a given grid 
the first case requires about 1100 iterations to converge, while 
subsequent cases for the same grid require about 300 iterations 
starting with the solution for the previous case. 

5 Results and Discussion 

Three experimental cases for the C3X vane, three for the 
VKI rotor, and two for the ACE rotor were analyzed for compar- 
ison. The values of various parameters for these cases are given 
in Table 1. In this table, the derived film-cooling parameters 
are based upon the assumption of one-dimensional compressible 
flow through the hole. For the C3X vane, cases 44135 and 
44155 represent the minimum and maximum blowing ratio, 
respectively, for the showerhead holes, while case 44355 repre- 
sents the wannest coolant. For the VKI rotor, case number is 
designated by the ratio of coolant mass to the main-flow mass; 
the main difference between the 3.09 and 3.32 percent cases 
being the colder coolant for the latter. For the ACE rotor, case 
6115 represents about half the blowing ratio for the case 6109. 
In all cases except case 6115 for the ACE rotor, the blowing 
ratio is high ( B r s 1.0). 


Figure 5 shows the coolant velocity and temperature profiles 
at the hole exit, the effect of which was studied on the heat 
transfer coefficient at various blade surfaces. The profiles shown 
in Fig. 5 are at the centerline through the hole, with the abscissa 
covering the hole. Clearly, the ?th power law profile does not 
satisfy the zero gradient condition at the centerline. The polyno- 
mial profile shown is for V m * 2.449 = twice the maximum 
velocity for the ?th power law profile, and V { - 0. It was found, 
however, that results for the heat transfer coefficient correspond- 
ing to the polynomial profile did not change more than ~5 
percent for -0.4 ^ V t ^ 0.4. 

We now present results for the three blades separately. Figure 
6 shows the nonuniform experimentally determined temperature 
on the C3X vane surface for the three cases. These temperature 
values were specified as the boundary condition for the blade 
surface temperature in the code as well. In this and later figures, 
except Figs. 10-12, s represents the normalized distance along 
the pressure or suction surface of the blade. Besides the some- 
what erratic temperature variation over \s\ > 0.25, caused by 
the internal cooling holes (not shown in Fig. 1) in the active 
part of the blade, there is a sharp drop in temperature at each 
end of the insulated portion of the blade (\s\ « 0.25). 

Figures 7—9 provide the effect of coolant velocity and tem- 
perature distribution at the hole exit on the normalized heat 
transfer coefficient ( solid and dash curves) at the blade surface 
in comparison with experimental data (Hylton et al., 1988), 
denoted by □, for the cases 44135, 44155, and 44355, respec- 
tively, at a spanwise location (near midspan) where the experi- 



Fig. 7 Effect of coolant velocity and temperature profiles on the normal- Fig. 9 Effect of coolant velocity and temperature profiles on the normal- 
ized heat transfer coefficient at the C3X vane surface (case 44135) ized heat transfer coefficient at the C3X vane surface (case 44355) 
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Fig. 10 Effect of coolant velocity and temperature profiles on the nor- pjg_ 12 Effect of coolant velocity and temperature profiles on the nor- 
malized heat transfer coefficient at the VKI rotor surface (2.07 percent malized heat transfer coefficient at the VKI rotor surface {3.32 percent 

coolant case) coolant case) 


mental data were taken. The heat transfer coefficient values in 
these figures have been normalized with respect to an arbitrary 
value, h 0 = 1135.6 kW/m 2 -K, as per Hylton et al. (1988). 
There are no data given for about 25 percent of surface length 
on either side of the leading edge since this portion contained 
the plenum chambers for injection of the colder gas and was 
insulated from the rest of the blade in the experimental tests 
(cf. Fig. 1). The fluctuations in the data are due to the nonuni- 
form blade surface temperature in the experimental data. The 
nine short vertical lines near the center-bottom of these figures 
denote the location of film cooling rows. For all three cases 
the normalized heat transfer coefficient corresponding to the 
polynomial profiles of coolant velocity and temperature at the 
hole exit is about 50-60 percent higher than that corresponding 
to the |th power law over most of the suction surface. However, 
the difference over the pressure surface is only about ±10 per- 
cent for the case 44135 (Fig. 7), about ±20 percent for the 
case 44155 (Fig. 8), and about 35 percent for the case 44355 
(Fig. 9). Also, the |th power law profile results seem to match 
better with the experimental data than the polynomial profile 
results for all the cases. 

For the C3X vane, we may point out that experimentally, the 
heat flux on the blade (inner) surface was calculated from a 
finite element analysis of conduction within the blade, and the 
heat transfer coefficient was then found by dividing this heat 
flux by (T 0 - T„). The present study computes the heat transfer 



Fig. 11 Effect of coolant velocity and temperature profiles on the nor- 
malized heat transfer coefficient at the VKI rotor surface (3.09 percent 
coolant case) 


coefficients on the blade (outer) surface from the three-dimen- 
sional Navier-Stokes analysis. Also, the uncertainty in the heat 
transfer coefficient measurement varies from about 10 percent 
near |s| » 0.3 to about 22 percent near \s\ ** 0.95 (Hylton et 
al., 1988). 

Let us now turn to the comparison of results for the VKI 
rotor. For this rotor, the blade was specified to be isothermal 
with TJT 0 = 0.7, Figures 10-12 display the effect of coolant 
velocity and temperature distribution at the hole exit on the 
normalized heat transfer coefficient (solid and dash curves) at 
the blade surface in comparison with experimental data ( Camci 
and Arts, 1990), denoted by □, for the three cases at midspan. 
In these figures, h n corresponds to the heat transfer coefficient 
at an uncooled blade surface. Also, the abscissa represents the 
surface distance along the blade normalized by the true chord, 
in conformity with the available experimental data (Camci and 
Arts, 1990). The six short vertical lines in the center-bottom 
of these figures denote the location of film-cooling rows. From 
these figures we find, in contrast to the results for the C3X 
vane, that the differences in h/h n values corresponding to the 
two coolant profiles at the hole exit are negligible over most of 
the suction surface and considerable (about 30 percent) over 
the pressure surface. The heat transfer coefficient corresponding 
to the polynomial profiles of coolant injection is generally larger 
than that corresponding to the ?th power law profiles. Also, 
comparison with experimental data is fair for either of the cool- 
ant profiles results. The experimental uncertainty in the mea- 
surement of heat transfer coefficient is about 5 percent except 
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Fig. 13 Blade surface temperature for the two cases (ACE rotor) 


Journal of Turbomachinery 


APRIL 1997, Vo!. 119/349 






Pressure Suction 


Fig. 14 Effect of coolant velocity and temperature profiles on the heat 
transfer coefficient at the ACE rotor surface (case 6115) 


in the shower-head region where it is estimated to be as high 
as 10- 15 percent (Camci and Arts, 1990). Due to normalization 
with respect to h ny however, each experimental data point repre- 
sents two measurements, one for the cooled and another for the 
uncooled blade. The same is true for the computations as well, 
and may be the reason for the relatively good comparison with 
the experimental data for the VKI rotor in contrast to that for 
the C3X vane in Figs. 7-9. We may also point out that the 
suction and pressure surface holes on the VKI rotor are conical 
while the showerhead holes on the VKI rotor as well as all 
holes on the C3X vane are cylindrical. This may account for 
the relatively smaller effect of coolant velocity and temperature 
distribution at the hole exit on the heat transfer coefficient at 
the VKI rotor surface in contrast to that for the C3X vane. 

We now turn to the discussion of results for the ACE rotor. 
Figure 13 shows the nonuniform experimentally determined 
temperature on the ACE rotor surface for the two cases ana- 
lyzed. These temperature values were specified as the boundary 
condition for the blade surface temperature in the code as well. 
Figures 14 and 15 show the effect of coolant velocity and tem- 
perature distribution at the hole exit on the heat transfer coeffi- 
cient (solid and dash curves) at the blade surface in comparison 
with experimental data (Norton et al., 1990), denoted by ■, 
for the two cases at midspan. The short vertical lines in these 
figures as well as Fig. 16 denote the location of film cooling 
rows. Due to the absence of showerhead cooling, this is the 
only case where Mayle’s transition criterion and Forest’s model 
for augmentation of laminar heat transfer at the leading edge 
due to free-stream turbulence are useful. Clearly, both models 
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Fig. 15 Effect of coolant velocity and temperature profiles on the heat 
transfer coefficient at the ACE rotor surface (case 6109) 



Fig. 16 Effect of Crawford’s model on the heat transfer coefficient at 
the ACE rotor surface (case 6109) 


seem to be effective in yielding results that compare well with 
the experimental data. While the effect of coolant profiles at 
the hole exit is negligible for the low blowing ratio case 6115 
in Fig. 14, it is considerable (as much as 50 percent higher h 
corresponding to the polynomial profile) over the pressure sur- 
face for the high blowing ratio case 6109 in Fig. 15. For both 
cases, however, the effect on the suction surface is negligible, 
in conformity with the results for the VKI rotor but in direct 
contrast to those for the C3X vane. The polynomial profile for 
the coolant velocity and temperature distribution at the hole exit 
results in a very good comparison with the experimental data 
downstream of the cooling holes on the pressure surface in Fig. 
15. Norton et al. (1990) could not get such a good comparison 
on the pressure surface; their prediction is similar to ours for 
the jth power law profile. 

Figure 16 shows the effect of Crawford’s model on the heat 
transfer coefficient at the ACE rotor surface for the case 6109. 
The results in Fig. 16 correspond to the polynomial profiles for 
coolant velocity and temperature at the hole exit. Clearly, the 
effect of Crawford’s model is felt only within the region of 
holes; it is negligible downstream of the holes. We may note that 
due to the absence of showerhead cooling, the flow upstream of 
the holes is essentially two dimensional (no spanwise variation) 
for the ACE rotor. 

6 Conclusions 

A three-dimensional Navier- Stokes code has been used to 
study the effect of coolant velocity and temperature distribution 
at the hole exit on the heat transfer coefficient on three film- 
cooled turbine blades, namely, the C3X vane, the VKI rotor, 
and the ACE rotor. Results are also compared with the experi- 
mental data for all the blades. Moreover, Mayle’s transition 
criterion (Mayle, 1991), Forest’s model for augmentation of 
leading edge heat transfer due to free-stream turbulence (Forest, 
1977), and Crawford’s model for augmentation of eddy viscos- 
ity due to film cooling (Crawford et al., 1980) are used. It is 
found that different velocity and temperature distributions of 
coolant at the hole exit can lead to as much as a 60 percent 
change in the heat transfer coefficient at the blade surface in 
some cases. Also, different effects are observed on the pressure 
and suction surface depending upon the blade as well as upon 
the hole shape, conical or cylindrical. Thus specification of 
proper conditions at the hole exit is important in film-cooling 
applications. This calls for a detailed analysis of the in-hole and 
near-hole physics under conditions relevant to the gas turbine. 
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